Interference-induced splitting of resonances in spontaneous emission 
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We study the resonance fluorescence from a coherently driven four-level atom in the Y-type con- 
figuration. The effects of quantum interference induced by spontaneous emission on the fluorescence 
properties of the atom are investigated. It is found that the quantum interference resulting from 
cascade emission decays of the atom leads to a splitting of resonances in the excited level populations 
calculated as a function of light detuning. For some parameters, interference assisted enhancement 
of inner sidebands and narrowing of central peaks may also occur in the fluorescence spectrum. We 
present a physical understanding of our numerical results using the dressed state description of the 
atom-light interaction. 

PACS numbers: 42.50.Ct,42.25.Hz,32.50.-f d 



I. INTRODUCTION 

The study of quantum interference effects in the spon- 
taneous emission of excited atoms has attracted substan- 
tial attention in the literature P, i, i, i, i, i, 0, B i, US] ■ 
The interference in spontaneous emission occurs when a 
pair of excited levels of an atom are coupled by the same 
vacuum modes to other levels. Many remarkable features 
have been predicted employing the mechanism of interfer- 
ences in the spontaneous emission of atoms [H, 0, H, I3j @ ■ 
The early work of Agarwal on this subject demonstrated 
trapping of populations in the degenerate excited levels 
of a V-type atom For a non-degenerate V system 
in free space, Zhu et al. predicted the existencjj of a 
dark line in the spontaneous emission spectrum [2]. By 
considering an open V system where the excited atomic 
levels are coupled by a coherent field to another auxiliary 
level, Scully, Zhu, and coworkers showed the possibility of 
spectral line elimination and spontaneous emission can- 
cellation Q via quantum interference. Phase dependent 
spectral narrowing Q and pulse propagation dynamics 
[5| have also been investigated using the four- level atomic 
model of Ref. @]. 

Since the fluorescence properties of a driven atomic 
system results from its spontaneous emission, studying 
the influence of interference in such processes has become 
an important topic of research 0,11, [1, [3 • The driven 
V system has been shown to exhibit many interference 
effects such as fluorescence quenching [6J, ultranarrow 
spectral lines [7] , anticorrelated photon emissions Q , en- 
hanced squeezing in the fluorescence field Q , and collec- 
tive population trapping [l^l . All these effects assume an 
existence of non-orthogonal dipole moments of the atomic 
transitions for the interference to occur [ij. However, in 
real atomic systems, it is difficult to meet this condition. 
Different schemes involving cavities with preselected po- 
larization [ll| ^oherent- and dc- field induced splitting of 
atomic levels [l^, [l3 bave been proposed later to bypass 
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the condition of non-orthogonal dipole moments. Fur- 
ther, the work on spontaneously generated interferences 
has been extended to four level atoms in different config- 
urations. The resonance fluorescence spectrum of driven 
four level atoms in the A- and V- type configurations has 
been extensively studied by Li et al. [13, [3] ■ Recently, 
Anton et al. [16| have examined a driven four level atom 
with three excited states and showed that a high popula- 
tion inversion may be achieved in the system due to the 
interference in spontaneous decay channels. 

In this paper, we consider a four-level atom in the Y- 
type scheme (as shown in Fig. 1) which was proposed ear- 
lier for studies on two photon absorption [13, ^M, ■ It is as- 
sumed that the excited atomic states are near-degenerate 
and decay spontaneously via the same vacuum modes to 
the intermediate state. The atom in the intermediate 
state can further decay to the ground state. Since the 
cascade decays (|1) ^ |3) -> |4) and |2) ^ |3) |4)) 
of the atom to its ground state from the two initially 
populated excited states lead to an emission of the same 
pair of photons, the quantum interference exists in decay 
processes. We investigate the role of the interference in 
the resonance fluorescence from the atom when driven by 
two coherent fields. 

The paper is arranged as follows. In Sec. II, we present 
the atomic density matrix equations, describing the in- 
teraction of a Y-type atom with two coherent fields, when 
the presence of quantum interference in decay channels is 
included. The population dynamics of the driven atom 
in the steady state is then studied in Sec. III. In Sec. 
IV, we analyze the fluorescence spectrum of the atom 
and identify the origin of interference effects using the 
dressed states of the atomic system. Finally, the main 
results are summarized in Sec. V. 



II. DRIVEN Y-TYPE ATOMIC SYSTEM AND 
ITS DENSITY MATRIX EQUATIONS 

We consider (Fig.l) the four-level Y-type atom having 
two closely lying excited states with the energy separa- 
tion hWi2. In this scheme, the excited atomic states 11) 
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where Aa = W13 — uja denotes the detuning between the 
atomic frequency (1^13) of the |1) |3) transition and 
the frequency of the apphed field Ea- Similarly, Af, = 
LO34 — ujb corresponds to the detuning of the field applied 
on the lower transition. 

We use the master equation framework to include re- 
laxation processes. With the inclusion of decay terms, 
the time evolution of the atomic density matrix describ- 
ing the atom-field interaction obeys 
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Pii 



-271P11 -I- i^i{p3i - P13) -pV7i72(/0i2 + P21), (5) 



FIG. 1: The level scheme of the Y-type atom driven by co- 
herent fields. 



P22 = -272/522 +if^2(p32 - P23) -PV7172(/012 +P21), (6) 



and 1 2) decay spontaneously to the intermediate state 
1 3) with rates 271 and 272, respectively. In addition, the 
atom in the intermediate state |3) can undergo sponta- 
neous emissions to the ground state |4) with decay rate 
273. We assume that direct transitions between the ex- 
cited states |1) 1 2) and that between the excited and 
ground states |2) |4) of the atom are forbidden in 
the dipole approximation. A coherent field of frequency 
uja (amplitude Ea)is set to couple the upper transitions 
|1),|2) ^ |3) and another field of frequency Wf, (ampli- 
tude Et) drives the lower transition |3) ^ |4). It is fur- 
ther assumed that the transition frequencies (0713,(^23) 
of the upper transitions arc widely different from that of 
the lower transition (lu^^). The Rabi frequencies of the 
atom-field interaction are represented as Qi — pLi^.Ea/h, 
^2 = ^23 -Ea/h, and ^3 = fl34.Eb/h with fimn being the 
dipole moment of the atomic transition from |m) to \n). 
The Hamiltonian of the atom-field interaction is given in 
the dipole and rotating wave approximations to be 



H = fvjJuAii -t- hjJ2/LA22 + hhJ34:A33 



-h{niAi3e-'^-* + n2A23e-"^-' + H.c.) 
-?i(n3^34e-'"''*-t-H.c.). 



(1) 



Here, the zero of energy is defined at the ground state 
|4), and hiUmn is the energy difference between the states 
|m) and \n). The operators Amn = |'tt-)(?t-| represent the 
atomic population operators for m — n and transition 
operators for m ^ n. The state \ip{t)) of the atomic 
system obeys the Schrodinger equation 



(2) 



It is helpful to use the interaction picture by making an 
unitary transformation \(j)) = eKp{iHot/h)\'ip) with 

Ho = h{bJa + LUb)Aii + h{uJa + UJh)A22 + hjJbA33. (3) 

In the interaction picture, the Schrodinger equation for 
the state \(t>) will have the effective Hamiltonian given by 

Hi = ;i(A„ + Afc)ylii + ;i(A„ + Ab - 14^12)^22 

-|-ftAfcyl33 - h{niAi3 + 172^23 + H.c.) 
-ft(r!3^34 + H.c.), (4) 



P33 = 271P11 -I- 272P22 - 273P33 -I- i^i{pi3 - P3i) 
+i^2{p23 - P32) + i^3ipi3 - P34) 

+2pV7i 72(^12 + P21), (7) 

P12 = -(71 + 72 + iWi2)pi2 + i^lP32 - i^2Pl3 

-P\/7i72(pii +^22), (8) 

PU = -(71 + 73 + «Aa)pi3 + i^l{p33 " Pll) " i^2Pl2 

~i^l3Pi4 -pV7i72 P23, (9) 

P23 = -[72 + 73 + «(Aa - Wl2)]p23 + «f^2(p33 " P22) 
-ifllP2l - ifl3P24 - PV7172 Pl3, (10) 

P34 = "(73 + «Ab)p34 -I- i03(p44 - P33) + iQipu 

+i^2P2i, (11) 



Pl4 = -[71 + K^a + A6)]pi4 + iQip3i - i^3Pl3 

-PV7i72 P24, (12) 

P24 = -[72 + «(Aq + Ab - H^12)]P24 + i^2P34 " if^3P23 

-pV7i72 Pi4, (13) 

In writing Eqs. ((5l)- (fT3| . we have assumed that the trace 
condition pn + P22 + P33 + P44 = 1 is obeyed. The cross- 
coupling term p = /Ii3./l23/|/Ii3||At23| arises due to the 
quantum interference in spontaneous decay transitions. 
This comes because the decays from the excited states 
|1) and 1 2) are coupled by the vacuum field. When p = 
±1, the interference effects are maximal, whereas if the 
dipoles are orthogonal (p = 0) there is no interference 
effect in spontaneous emission. 

The density matrix equations ((5|l- (|13p can be rewritten 
in a more compact matrix-form by the definition 
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* = (Pll , P22 , P33 , /O12 , Pis , P23. , Pl4 , P24 , P34, P2I , P31 , P32 , p41 , P42 , P43) 



(14) 



r 



Substituting Eq. HI]) into Eqs. ©-([ni), we get the 
matrix equation for the variables {t) 



— * = + 
at 



(15) 



where j is the j-th component of the column vector ^ 
and the inhomogeneous term / is also a column vector 
with non-zero components 



/i5 -zfia 



(16) 



In Eq. (1151) . ^ 15x15 matrix whose elements are 
time independent and can be found explicitly from Eqs. 
([5))- P^ . The steady state solutions of the density ma- 
trix elements can be found by setting the time derivative 
equal to zero in Eq. (|15p : 



(17) 



III. STEADY STATE POPULATIONS 

We first study the population dynamics of the driven 
atom in steady state using Eq. (fTT]) . In Fig. 2, 
we show the excited and intermediate level populations 
[pii(oo) = pii,p22(oo) = P22>P33(oo) = P33] vcrsus the 
detuning Aq for different decay rates. All the frequency 
parameters such as decay rates, detuning, and Rabi fre- 
quencies are scaled in units of 73. It can be seen in Fig. 
2 that interference effects [p — 1) are less prominent for 
7i,72 < 73. This feature is expected as the interference 
terms scale as p^7i72 in Eqs. ([5l)- (fT3l) . Further, the 
graphs show that the excited level populations exhibit a 
resonance at the value of detuning close to « in 
the absence of interference (p = 0). More generally, the 
resonances in excited level populations and P22 occur 
when the two photon resonance conditions A^ -f A;, = 
and Aa -I- Af, = 1^12 for the |1) ^ |4) and |2) ^ |4) 
transitions are respectively satisfied [19|. The effect of 
interference is seen to enhance little the population in the 
excited atomic state when the one photon transitions are 
resonant, A^ = 0, A^ = [see Fig. 2(a)]. Interestingly, 
for the case of 71,72 ^ 73, the interference leads to a 
splitting of resonances in the excited level populations as 
shown in Fig. 2(b). This result is purely the effect of 
couplings between the different decay pathways that the 
excited atom can take. It should be borne in mind that 
both one- (pis, P23) and two-photon (/9i4, P24) coherences 
contribute in the interference among the decay pathways. 



To explore further the interference induced splittings 
of population resonances, we consider the case of near- 
degenerate excited levels (M^i2 ^ 71,72) and take the 
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FIG. 2: Steady state population of atomic levels as a function 
of the detuning Aa for the parameters 73 = 1, W12 ~ 5, 
At = 0, ^1 = f72 = fis = 3, 71 = 72 = 0.5 (a), 71 = 72 = 2 
(b). Actual values of P33 are three [six] times than shown in 
(a) [(b)]. The solid (dashed) curves are for p — 1 {p = 0). The 
curves for (not shown) have a similar behavior as that of 
Pii- 



high intensity limit {fli,fl2,fl3 > W^i2, 7i, 72, 73) of ap- 
plied lasers. For simplicity, we assume equal decay rates 
7i = 72 = 7 for the upper transitions and examine two 
different cases, (a) 7 ^ 73, (b) 7 = 73, with respect to 
the decay rate of the lower transition. The numerical re- 
sults are shown in Fig. 3 which are to be compared with 
Fig. 2. It is found that resonances in excited level pop- 
ulations occur at Aq = ztil [see Fig. 3(a)] in the limit 
^ 7 S> 73, where f^i = f72 = f^3 = is considered. In 
the case of equal decay rates 7 = 73, analytical expres- 
sions for the population (pn) can be obtained compactly 
in the presence {p = 1) and absence {p = 0) of interfer- 
ence as 
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_ 41^4 Ag + 4»«A4, + 40171" Ag + 160»i" 

^"''^ ^ 8172 A8 + 81^6 A6 + 64178 A4 + 2401^10 A2 + 9601^10 ' >^ «^ 



where all the parameters have been scaled in units of 7. 
These analytical formulas account well for the numerical 
results in Fig. 3(b). In order to understand physically 
the effect of interference, the atomic dynamics is further 
studied in the bases of symmetric and anti-symmetric 
states Q defined by 



1 



(V7^|l>+^/^|2)), 
(V7^|l)-V7r|2)). (19) 



V71 + 72 
V7i + 72 

With 7i = 72 = 7 and using Eq. p9)) . the Hamiltonian 
Eq. (|3|) can be rewritten as 

Hr= fi(A,--^)(|.)(.| + |a)(a|) 

+ ^(k>(«l + - nV2ni\s){3\ + \3){s\) 



-?in(|3)(4| + |4)(3|). 



(20) 




From the above Hamiltonian, it is seen that only the sym- 
metric state \s) is interacting with the light field. How- 
ever, the antisymmetric state \a) may be populated by 
its coupling with the symmetric state because of the sep- 
aration energy (hWu) between the excited atomic levels. 
This can become clear by analyzing the density matrix 
equations in the bases (fTT 



Paa 
Pss 



iWi2 _ iWi2 

~Pas Psa-! 



2 

-47P, 



2 

iWi2 



-Pas 



2 2 
iV2n{p3s - Psz)- 



(21) 



Here, the case of maximal quantum interference {p — 1) 
has been assumed. It is evident that the antisymmetric 
state is a non-decaying state for p = \ and it is coupled 
to the symmetric state for W12 7^ (though small as in 
Fig. 3). In Fig. 4, the steady state populations of the 
symmetric and antisymmetric states are plotted for the 
same parameters of Fig. 3(a). The graphs show that the 
splitting of resonances occurs due to a high population 
of the antisymmetric state. We have so far assumed a 
fixed value for the lower transition detuning (A^ = 0) 
and studied the dependence of populations on the up- 
per transition detuning (A^). However, the results (not 
shown) will be qualitatively similar even in the general 
case of varying both A^ and A^. 
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FIG. 3: Steady state population of atomic levels as a function 
of the detuning Aa for the parameters 73 = 1, VK12 = 0.2, 
Ab = 0, 17i = 172 = = 10, 71 = 72 = 5 (a), 71 = 72 = 1 
(b). Actual values of P33 are six times than shown. The solid 
(dashed) curves are for p = 1 (p = 0) . 



FIG. 4: Steady state populations p^^ and p^^ as a function of 
the detuning Aa for the same parameters of Fig. 3(a) with 
0=1. 
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FIG. 5: The incoherent spectrum of the fluorescent field gen- 
erated by the dipoles (a) pL^' and (b) pL^'' for the parameters 

73 = 1, Wl2 = 10, Aa = Ab = 0, ^1 = ^2 = 10, ^3 = 5, 

7i = 72 = 3. The solid (dashed) curves are for p = 1 (p = 0). 



IV. RESONANCE FLUORESCENCE 
SPECTRUM 



We now proceed to the study of the resonance fluo- 
rescence spectra of the driven atom. Since the atom is 
driven by two coherent fields, each field induces its own 
atomic dipole moment which then generates a scattered 
field. However, the fields scattered by the upper- and 
lower- transitions in the atom will have no correlations 
because the applied fields [Ea ,^^&) arc of quite different 
carrier frequencies (wa, Wf,). In the interaction picture, 
the negative- and positive-frequency parts of the polar- 
ization operators are written as 



Mi3e*""*|l>(3|+/J23e*""*|2)(3|, 



(22) 



or equivalently the real part of its Laplace transform: 

/•OO 

Sa{uj) = Re / hm (P£;)(t + r).P£+) (0)e-'"^dr, 

f'OO 

5,(c^)=Re / lim(p(;)(t + r).PLt)(i))e-"^dr.(24) 

Jo t-»oo 

Here, the index a (b) refers to the spectrum of the flu- 
orescence light emitted by the atom with a central fre- 
quency LOa i^b)- The Laplace transformation with vari- 
able Z — iuj oi the correlation function, defined in the 
spectrum above, has a pole at Z = iua iZ — iuji,) which 
attributes to the coherent Rayleigh scattering of the spec- 
trum. The incoherent part is obtained by removing the 
contributions of the poles. 

With the application of the quantum regression the- 
orem EE and using the steady state solutions 
Eq. P7|) of the density matrix elements, the incoherent 
fluorescence spectra can be obtained as 

Sa{0j) = Rc||/Ji3|^[Mll,9Pi4+Mii^3Pi3-HMii42Pl2 
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+M1IMP11 + NiijIjPis] + /J23-/Jt3 

X [^12.9Pl4 + Mi2,3Pl3 + ^^12,12Pl2 + A^12,ll 
15 

XPll + ^12j-^iPl3] + M13-M23 [^^11,11P21 

J = l 

+Mll^9P2i + -^Al.3P23 + ^'^11.12P22 (25) 
15 

+ X!^11J^J^23] + lM23|^[Afl2,llP21 +^^2,9 
J = l 

15 

XP24 + ^^12,3P23 + A^12,12P22 + Y ^^"^^oh 



'P23. 



where the matrices M = {Z — L) \z=i(uj-uja) ^^-^ N = 
L^^M. Similarly, 

Sb(uj) = Re||/i34|^ [Mi5,13P31 + ^'^15,14]032 + ^^15,15^33 



10 



(26) 



PL;^W = /Z34e^"^*|3)(4|, pW(t) = [PL;^^)]^ (23) 

To calculate the fluorescence spectra, we need the two- 
time expectation values of the polarization operator. 
The spectrum of resonance fluorescence is defined by 
the Fourier transformation of the two-time correlation 



with the matrices M = {Z — L) \z=i(Lj~i^b) and N — 

The set of equations (P5|) and can be used to ob- 
tain numerically the spectral characteristics of the driven 
atom. Figure 5 displays the numerical results by assum- 
ing equal decay rates 71 — 72 of the upper atomic tran- 
sitions. The spectra Sa{oj) and S'h(w) are scaled in units 
of |/ii3p73'^ ^'^^ lA'34p73'^ respectively. In the presence 
of quantum interference {p = 1), the spectrum shows 
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FIG. 6: The dressed state eigenvalues versus the Rabi fre- 
quency f2 for the parameters 73 — 1, W12 = 10, Aa = At, = 0, 
^3 = 5. 



the typical line narrowing effect, as discussed in earlier 
publications 0, H, , in the fluorescent field with the 
central frequency uJa [see Fig. 5(a)]. However, the spec- 
tral features get remarkably modified in the fluorescent 
field emitted by the lower atomic transitions. It is seen 
that the inner sideband in the fluorescence spectrum gets 
enhanced due to interference with a corresponding re- 
duction in the intensity of the outer sidebands [compare 
solid and dashed curves in Fig. 5(b)]. A physical un- 
derstanding of this interesting result can be obtained in 
the dressed state description of the atom-field interac- 
tion. The dressed states are defined as eigenstates of the 
time independent Hamiltonian ^ : 



Hj\<S>) = hX\<i>) 



(27) 



In the general parametric conditions, it is difficult to find 
analytical solutions to the eigenvalue equation (|27p . For 
simplicity, the case of two photon resonance — ~ 
is assumed in the following. In this case, there exists an 
eigenstate \d) with the eigenvalue = as 



\d) = 



1 



[r!3|l)-r!i|4)]. 



(28) 



The non-zero eigenvalues and the corresponding eigen- 
states can be obtained by diagonalizing the Hamilto- 
nian Hj in the basis of bare atomic states. We con- 
sider a special choice of parameters ili = fl2 = and 
^3 = 1^12/2 [as in Fig. 5] which allows for simple 
analytical solutions. The operator Hj has eigenstates 
|m), |±) with eigenvalues (in units of h) Am = — M^i2/2, 

A± = (-1^12 ± 
where 



W?2 



32(f72 



4 



))/4, respectively. 



1 



n\i)-n\2) + ^\3) 



l±> = ^± 



m) 



2 



|4) 



(29) 



{W12 + A±) 



14^19 

|2)-A±|3) + ^|4) 




FIG. 7: Steady state population of dressed states, (a) pmm, 
pdd and (b) p++, p — , as a function of the Rabi frequency Q 
for the parameters 73 — 1, W12 = 10, Aa = Aj = 0, f23 = 5, 
71 = 72 = 3. The solid (dashed) curves are for p — 1 {p — 0). 



with A^± = l/^A^[l + Q^/{Wi2 + A±)1 + + 

In order to interpret the numerical results in Fig. 5, we 
study the behavior of the dressed states in steady state 
with the inclusion of decay processes using Eq. ([TT]). In 
Figs. 6 and 7, the dressed state eigenvalues (Avalues) and 
its populations are shown as a function of the Rabi fre- 
quency for the fixed value of Sla = W12/2. Note that 
the eigenvalues Xd and Am are independent of the pa- 
rameter ri. The peaks in the fluorescence spectrum can 
be attributed to transitions between the dressed states 
1$') ^ 1$) ($,$' ^d,m,+,-). Forp = 0andf} >M^i2, 
the dressed states Id), |-|-), and |— ) are well populated as 
shown in Fig. 7. The fluorescence peaks in Fig. 5 occur 
at the energy differences between these states. However, 
in the presence of interference {p = 1), the atomic pop- 
ulation accumulates mostly in the dressed state |m) [see 
Fig. 7(a)]. This can be explained as due to a destruc- 
tive quantum interference among the spontaneous decay 
pathways. The rate of transitions |<I>) |<I>') between 
dressed states |$) and |$') is given by the squared dipole 
matrix elements which for the emission lines with the 
central frequencies uJa and LUb becomes 



R 



where 



|M13| 



l(1>'|P£t^l*)l 

\tiiz?[Cl^ + C|* + 2pCi$C2*]C2$,, (30) 

i(i>'ip£-:^i'i>>i' 

|/234pC'3$C|$,, (31) 

|/i23| has been assumed and Ci$'s [i = 



7 



1, 2, 3, 4] denote the coefficients of the bare atomic states 
\i) in the dressed state 1$). As seen in the above Eq. 
pop , there is an interference term (p term) in the tran- 
sition rate i?| ^, as a result of spontaneous decays along 
the |1) |3) and |2) |3) transitions. For the dressed 
state |m), this term cancels the square factors in the case 
of maximal interference (p = 1), thus suppressing the 
atomic decay, narrowing the spectral lines [shown in Fig. 
5(a)] and enhancing the population [shown in Fig. 7(a)] 
in this state. However, the dressed state |m) can decay 
because of spontaneous emissions along the |3) |4) 
transitions even when p = 1 [see Eq. (|3ip ] . This leads to 
the enhancement of the inner sideband in the spectrum 
[shown in Fig. 5(b)] of the fluorescence light emitted 
by the lower transitions in the atom. It is because only 
the states |m) and \d) are populated mainly in steady 
state. Finally, we note that the existence of atomic steady 
state and discussions so far assume the non-degenerate 
(Wi2 7^ 0) case of excited atomic levels. In the degen- 
erate case {Wi2 — 0), there exists no unique solution 
to Eq. I|15p in steady state. In fact, the steady state 
fluorescence properties become dependent on the initial 
conditions due to degeneracy of the dressed states of the 
Hamiltonian. 



V. SUMMARY 

We have investigated the resonance fluorescence from 
a driven Y-type atom when the presence of interference 
in spontaneous decay channels is important. At first, the 
steady state dynamics of the atom was studied using the 
density matrix approach. We have shown that the decay- 
induced interference can lead to splitting of resonances in 
the excited level populations calculated as a function of 
light detuning. This has been explained as due to high 
population of a non-decaying anti-symmetric state of the 
atom. Then, the role of interference in the spectral char- 
acteristics of the driven atom was examined. It is found 
that the interference results in narrowing of central peaks 
and enhancement of inner sidebands in the fluorescence 
spectrum. A physical understanding of the numerical re- 
sults has been presented based on the dressed state the- 
ory of atom-field interaction. Clearly, the present work 
is open ended with the effects of interference in driven Y 
systems on two photon correlations and squeezing spec- 
tra remaining unexplored. Detailed investigations of such 
studies will be published elsewhere. 
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